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SYSTEM FOR LIMITING LOUDSPEAKER DISPLACEMENT 



Field of the Invention 

This invention generally relates to electro-acoustical transducers 
(loudspeakers), and more specifically to signal processing for limiting a vibration 
5 displacement of a coil-diaphragm assembly in said loudspeakers. 

Background of the Invention 

The Problem formulation 

A signal driving a loudspeaker must remain below a certain limit. If the signal 
is too high, the loudspeaker will generate nonlinear distortions or will be irreparably 
10 damaged. One cause of this nonlinear distortion or damage is an excess vibration 

displacement of a diaphragm-coil assembly of the loudspeaker. To prevent nonlinear 
distortion or damage, this displacement must be limited. 

Displacement limiting can be implemented by continuously monitoring the 
displacement by a suitable vibration sensor, and attenuating the input signal if the 
15 monitored displacement is larger than the known safe limit. This approach is 

generally unpractical due to the expensive equipment required for measuring the 
vibration displacement. Thus some type of a predictive, model-based approach is 
needed. 

Prior art solutions 

20 The prior art of the displacement limiting can be put into three categories: 



1. 



Variable cut-off frequency filters driven by displacement predictors. 



2. 



Feedback loop attenuators. 



3. 



Multi-frequency band dynamic range controllers. 



The prior art in the first category has the longest history. The first such system 
was disclosed in US Patent No. 4,1 13,983, "Input Filtering Apparatus for 
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Loudspeakers", by P. F. Steel. Further refinements were disclosed in US Patent No. 
4,327,250, "Dynamic Speaker Equalizer", by D. R. von Recklinghausen and in US 
Patent No. 5,481,617, "Loudspeaker Arrangement with Frequency Dependent 
Amplitude Regulations" by E. Bjerre. The essence of the prior art in the first category, 
5 utilizing a variable high pass filter with a feedback control for said displacement 
limiting, is shown in Figure la. 

In this category of loudspeaker protection systems (as shown in Figure la), a 
high-pass filter 12 of a signal processor 10 filters the input electro-acoustical signal 
22. Then a filtered output signal 24 of said high-pass filter 12 is sent to a loudspeaker 

10 20 (typically, through a power amplifier 18) and also fed to a feedback displacement 
predictor block 14. If the value of the displacement exceeds some predefined 
threshold value, a feedback displacement prediction signal 26 from the block 14 
indicated that and a cut-off frequency of the high-pass filter 12 is increased based on 
the feedback frequency parameter signal 28 provided to the high-pass filter 12 by a 

15 feedback parameter calculator 16 in response to said feedback displacement 

prediction signal 26. By increasing the cut-off frequency of the high-pass filterl2, 
lower frequencies in the input signal, which generally are the cause of the excess 
displacement, are attenuated, and the excess displacement is thereby prevented. 

The prior art in the first category has several difficulties. The high-pass filter 
20 12 and the feedback displacement predictor block 14 have finite reaction times; these 
finite reaction times prevent the displacement predictor block 14 from reacting with 
sufficient speed to fast transients. Bjerre presented a solution to this problem in US 
Patent No. 5,481,617 at the expense of significantly complicating the implementation 
of the displacement limiting system. An additional problem comes from the fact that 
25 the acoustic response of the loudspeaker naturally has a high-pass response 

characteristic: adding an additional high-pass filter in the signal chain in the signal 
processor 10 increases the order of the low-frequency roll-off. This can be corrected 
by adding to the signal processor a low-frequency boosting filter after the high-pass 
filter, as was disclosed by Steel in US Patent No. 4,1 13,983. However, this further 
30 complicates the implementation of the signal processing. 
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Prior art in the second category was disclosed in US Patent No. 5,577,126, 
"Overload Protection Circuit for Transducers", by W. Klippel. Figure lb shows the 
essence of a loudspeaker protection system describing this category. The output of 
the displacement predictor is fed-back into the input signal, according to a feedback 
5 parameter k, calculated by a threshold calculator. This category of the vibration 

displacement protection is simpler than the first category system described above, in 
that it does not require a separate high-pass filter. 

Prior art in the second category can be effective for the vibration displacement 
limiting. However, the feedback loop has an irregular behaviour around a threshold 

10 value, due to a modification of the loudspeaker's Q- factor, and an amplification at 
low frequencies. These effects can cause subjectively objectionable artifacts. In the 
above-mentioned US Patent No. 5,577,126, Klippel describes one solution to this 
problem: the attenuation of the signal processor is somewhat better behaved if the 
pure feedback signal path 16 is differentiated, as shown in Figure 3 of US Patent No. 

1 5 5,577,126. However, this causes significant and unnecessary attenuation of the higher 
frequency band. Therefore, signals that are not responsible for the excess 
displacement are likely to be attenuated, degrading the performance of the 
loudspeaker system. 

Prior art in the third category was disclosed in WO Patent Application No. 

20 PCT/EP00/05962 (International Publication Number WO 01/03466 A2), 

"Loudspeaker Protection System Having Frequency Band Selective Audio Power 
Control", by R. Aarts. Figure lc shows the essence of the third category loudspeaker 
protection system. The input signal is divided into N frequency bands by a bank of 
band-pass filters. The signal level in the n th frequency band is modified by a variable 

25 gain g n . The signals in the N frequency bands are summed together, and sent to the 
power amplifier and loudspeaker. An information processor monitors the signal level 
in each frequency band, as modified by each of the variable gains gi, g 2 , . . .g n . The 
information processor modifies the variable gains gi, g2 5 . . .g n in such a way as to 
prevent the excess displacement in the loudspeaker. The advantage of the third 

30 category approach is that the signal is attenuated in only that frequency band that is 
likely to cause the excess loudspeaker diaphragm-coil displacement. The remaining 
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frequency bands are unaffected, thereby minimizing the effects of the displacement 
limiting on the complete audio signal. 

The disadvantage of the third category displacement limiter is that there are no 
formal rules describing how the information processor should operate. Specifically, 
5 no formal methods are available for describing how the information processor should 
modify the gains g n so as to prevent the output signal from driving the loudspeaker's 
diaphragm-coil assembly to the excess displacement. The information processor can 
only be designed and tuned heuristically, i.e., by a trial-and-error. This generally 
leads to a long development time and an unpredictable performance. 

10 Summary of the Invention 

The object of the present invention is to provide a novel method of signal 
processing for limiting a vibration displacement of a coil-diaphragm assembly in 
electro-acoustical transducers (loudspeakers). 

According to a first aspect of the invention, a method for limiting a vibration 
15 displacement of an electro-acoustical transducer comprises the steps of: providing an 
input electro-acoustical signal to a low frequency shelving and notch filter and to a 
displacement predictor block; generating a displacement prediction signal by said 
displacement predictor block based on a predetermined criterion in response to said 
input electro-acoustical signal and providing said displacement prediction signal to a 
20 parameter calculator; and generating a parameter signal by said parameter calculator 
in response to said displacement prediction signal and providing said parameter signal 
to said low frequency shelving and notch filter for generating an output signal and 
further providing said output signal to said electro-acoustical transducer thus limiting 
said vibration displacement. 

25 According further to the first aspect of the invention, the electro-acoustical 

transducer may be a loudspeaker. 

Further according to the first aspect of the invention, the low frequency 
shelving and notch filter may be a second order filter with a z-domain transfer 
function given by 
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l + a u z +a 2t z 

wherein <r c is a characteristic sensitivity of the low frequency shelving and notch 
filter, bi. c and b 2 . c are feedforward coefficients defining target zero locations, and ai. t 
and a2- t are feedback coefficients defining target pole locations. Further, said 
5 parameter signal may include said characteristic sensitivity a c and said feedback 
coefficients a^and ai. t . 

Still further according to the first aspect of the invention, the method may 
further comprise the step of: generating said output signal by the low frequency 
shelving and notch filter. Further, the method may further comprise the step of: 
10 providing the output signal to said electro-acoustical transducer. Yet further, the 

output signal may be amplified using a power amplifier prior to providing said output 
signal to said electro-acoustical transducer. 

According further to the first aspect of the invention, the displacement 
prediction signal may be provided to a peak detector of the parameter calculator. Still 

15 further, after the step of generating the displacement prediction signal, the method 

may further comprise the step of: generating a peak displacement prediction signal by 
the peak detector and providing said peak displacement prediction signal to a 
shelving frequency calculator of the parameter calculator. Yet still further, the method 
may further comprise the step of: generating a shelving frequency signal by the 

20 shelving frequency calculator based on a predetermined criterion and providing said 
shelving frequency signal to a sensitivity and coefficient calculator of the parameter 
calculator for generating, based on said shelving frequency signal, the parameter 
signal. 

According still further to the first aspect of the invention, the input electro- 
25 acoustical signal may be a digital signal. 

According further still to the first aspect of the invention, said low frequency 
shelving and notch filter may be a second order filter with an s-domain transfer 
function given by 
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wherein Q c is a coefficient corresponding to a Q-factor of the electro-acoustical 
transducer, co c is a resonance frequency of the electro-acoustical transducer mounted 
in an enclosure, Q t is a coefficient corresponding to a target equalized Q-factor, co, is 
5 a target equalized cut-off frequency. Still further, Q c may be equal to 1/ 4l , when the 
electro-acoustical transducer is critically damped. Yet further, Q c may be a finite 
number larger than 1/ 4l , when the electro-acoustical transducer is under-damped. 

According to a second aspect of the invention, a computer program product 
comprising: a computer readable storage structure embodying computer program code 
10 thereon for execution by a computer processor with said computer program code, 

characterized in that it includes instructions for performing the steps of the first aspect 
of the invention indicated as being performed by the displacement predictor block or 
by the parameter calculator or by both the displacement predictor block and the 
parameter calculator. 

1 5 According to a third aspect of the invention, a signal processor for limiting a 

vibration displacement of an electro-acoustical transducer comprises: a low frequency 
shelving and notch filter, responsive to an input electro-acoustical signal and to a 
parameter signal, for providing an output signal to said loudspeaker thus limiting said 
vibration displacement of said electro-acoustical transducer; a displacement predictor 

20 block, responsive to said input electro-acoustical signal, for providing a displacement 
prediction signal; and a parameter calculator, responsive to said displacement 
prediction signal, for providing the parameter signal. 

According further to the third aspect of the invention, the parameter calculator 
block may comprise: a peak detector, responsive to the displacement prediction 
25 signal, for providing a peak displacement prediction signal; a shelving frequency 
calculator, responsive to the peak displacement prediction signal; for providing a 
shelving frequency signal; and a sensitivity and coefficient calculator, responsive to 
said shelving frequency signal, for providing the parameter signal. Further still, said 
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low frequency shelving and notch filter may be a second order digital filter with a z- 
domain transfer function given by 



-2 



l + a w z + # 2 -, z 

wherein cs c is a characteristic sensitivity of the low frequency shelving and notch 
5 filter, bi. c and b2. c are feedforward coefficients defining target zero locations, and ai. t 
and a 2 . t are feedback coefficients defining target pole locations. Yet further, said 
parameter signal may include said characteristic sensitivity <j c and said feedback 
coefficients ai . t and ai . t . 

Further according to the third aspect of the invention, the output signal may be 
10 provided to said electro-acoustical transducer or said the output signal is amplified 

using a power amplifier prior to providing said output signal to said electro-acoustical 
transducer. 

Still further according to the third aspect of the invention, the input electro- 
acoustical signal may be a digital signal. 

15 According further to the third aspect of the invention, the low frequency 

shelving and notch filter may be a second order filter with an s-domain transfer 
function given by 

H ,,_ s2 +sv c /Q c +to 2 c 
s* +s(d t /Q t +co; 

wherein Q c is a coefficient corresponding to a Q-factor of the electro-acoustical 
20 transducer, co c is a resonance frequency of the electro-acoustical transducer mounted 
in an enclosure, Q t is a coefficient corresponding to a target equalized Q-factor, o, is 
a target equalized cut-off frequency. Further, Q c may be equal to 1/ V2 , when the 
electro-acoustical transducer is critically damped. Yet still further, Q c may be a finite 
number larger than 1/ 4l , when the electro-acoustical transducer is under-damped. 

25 According still further to the third aspect of the invention, the electro- 

acoustical transducer may be a loudspeaker. 
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Brief Description of the Drawings 

For a better understanding of the nature and objects of the present invention, 
reference is made to the following detailed description taken in conjunction with the 
following drawings, in which: 

5 Figures la, lb and lc show examples of a signal processor and a loudspeaker 

arrangement for a first, second and third category signal processing systems for a 
loudspeaker protection (vibration displacement limiting), respectively, according to 
the prior art. 

Figures 2a shows an example of a signal processor with a loudspeaker 
10 arrangement utilizing a variable low-frequency shelving and notch filter driven by a 
feedforward control using a displacement predictor block, according to the present 
invention. 

Figures 2b shows an example of a parameter calculator used in the example of 
Figure 2a, according to the present invention. 

15 Figure 3 shows an example of response curves of a low- frequency shelving 

and notch filter (without a notch and Q c =0.707) for a critically damped loudspeaker, 
according to the present invention. 

Figure 4a and 4b show examples of displacement response curves for a 
loudspeaker which is critically damped and under-damped, respectively, by utilizing a 
20 low-frequency shelving and notch filter of Figure 3, according to the present 
invention. 

Figure 5a shows an example of response curves of a low-frequency shelving 
and notch filter (with a notch and Q c =6.4) for an under-damped loudspeaker, 
according to the present invention. 

25 Figure 5b shows an example of displacement response curves for a 

loudspeaker which is under-damped by utilizing a low-frequency shelving and notch 
filter of Figure 5a, according to the present invention. 

Figure 6 is a flow chart demonstrating a performance of a signal processor 
with a loudspeaker arrangement utilizing a variable low-frequency shelving and notch 
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filter driven by a feedforward control using a displacement predictor block, according 
to the present invention. 

Best Mode for Carrying Out the Invention 

The present invention provides a novel method for signal processing limiting 
and controlling a vibration displacement of a coil-diaphragm assembly in electro- 
acoustical transducers (loudspeakers). The electro-acoustical transducers are devices 
for converting an electrical or digital audio signal into an acoustical signal. For 
example, the invention relates specifically to a moving coil of the loudspeakers. 

The problems of the prior art methods described above for the displacement 
limiting is solved by starting with the first category approach, and making the 
following modifications: 

• Replacing the variable high-pass filter 12 (see Figure la) with a variable low- 
frequency shelving and notch (LFSN) filter; 

• Using a feedforward instead of a feedback control of the filter 12 by the 
displacement predictor block; 

• Employing a digital implementation; 

• Approximating the exact formulas for calculating required coefficients by 
finite polynomial series. 

According to the present invention, a signal processor with the above 
characteristics or a combination of some of these characteristics provides a 
straightforward and efficient system for said displacement limiting. Large signals that 
can drive the loudspeaker into an excess displacement are attenuated at low 
frequencies. Higher-frequency signals that do not overdrive the loudspeaker can be 
simultaneously reproduced unaffected. The behaviour of the limiting system can be 
known from its base operating parameters, and can therefore be tuned based on the 
known properties of the loudspeaker. 

Figure 2 shows one example among others of a signal processor with a 
loudspeaker arrangement utilizing a low- frequency shelving and notch (LFSN) filter 11 
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driven by a feedforward control using a displacement predictor block 14a for limiting 
a vibration displacement of an electro-acoustical transducer (loudspeaker) 20, 
according to the present invention. The limiting of the vibration displacement is 
achieved by modifying a transfer function of the LFSN filter 11 based on the output 
5 of the displacement predictor block 14a. 

As in Figure la, the LFSN filter 11 of a signal processor 10a filters the input 
electro-acoustical signal 22. Said input electro-acoustical signal 22 can be a digital 
signal, according to the present invention. Then a filtered output signal 24a of said 
high-pass filter 11 is sent to a loudspeaker 20 (typically, through a power amplifier 

10 18). But, according to the present invention, the input electro-acoustical signal 22 is 
also fed to a displacement predictor block 14a. If the value of the vibration 
displacement exceeds a predefined threshold value (that is a predetermined criterion), 
a displacement prediction signal 26a from the block 14a is generated and provided to 
the parameter calculator 16 which generates a parameter signal 28a in response to that 

15 signal 26a and then said parameter signal 28a is provided to the LFSN filter 11. Based 
on said parameter signal 28a, the transfer function of said LFSN filter 11 is modified 
appropriately and the output signal 24a of said LFSN filter 1 1 has the vibration 
displacement component attenuated based on said predetermined criterion. 

The LFSN filter 11 attenuates only low frequencies, which are the dominant 
20 sources of a large vibration displacement. The diaphragm-coil displacement can be 
predicted from the input signal 22 by the displacement predictor block 14a 
implemented as a digital filter. Generally, the required order of said digital filter is 
twice that of the number of mechanical degrees of freedom in the loudspeaker 20. The 
output of this filter is the instantaneous displacement of the diaphragm-coil assembly 
25 of the loudspeaker 20. The performance of the displacement predictor block 14a is 
known in the art and is, e.g., equivalent to the performance of the part 9 shown in 
Figure 2 of US Patent No. 4,327,250, "Dynamic Speaker Equalizer", by D. R. von 
Recklinghausen. Detailed description of the parameter calculator la is shown in an 
example of Figure 2b and discussed in detail later in the text. 
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The LFSN filter 11 can be designed, according to the present invention, as a 
second-order filter with an s-domain transfer fixnction given by 

^W=- ^ — r t 1 * 

5 + S(D t /Q t + CD, 

wherein £?c is a coefficient corresponding to a Q-factor (of the loudspeaker 20), co c is a 
5 resonance frequency of a loudspeaker 20 mounted in a cabinet (enclosure), in rad/s, Q t 
is a coefficient corresponding to a target equalized Q-factor, co, is a target equalized 
cut-off frequency (shelving frequency), in rad/s. The magnitude of the frequency 
response of the filter 11, a low-frequency gain, equals to co^/co 2 . Typical gain curves 

for this low-frequency shelving and notch filter 11 with Q c = Qt =1/V2 (the 
10 loudspeaker 20 is critically damped and the LFSN filter 11 does not have a notch) are 
shown in Figure 3 for five values of a>V w 2 c ratio. The ability of the LFSN filter 11 to 
limit the displacement is made clear in Figure 4a. 

Figure 4a shows an example among others of displacement response curves 
for the loudspeaker 20, which is critically damped by utilizing the LFSN filter 11 of 
15 Figure 3, according to the present invention. As the value of a>, is increased, the 

displacement response is attenuated as seen in Figure 4a. In the low frequency limit, 
the amount of attenuation varies as a> 2 . The mathematical detail behind this is 
discussed below. These displacement response curves are for a "critically damped" 
loudspeaker, i.e., one tuned to a Butterworth alignment (Q c = Q t = l/V2 ). 

20 Inexpensive loudspeakers often have an under-damped response, i.e., having 

values of Q c and Q t greater than l/ V2 . Figure 4b shows an example of displacement 
response curves for the loudspeaker 20 which is under-damped, by utilizing the LFSN 
filter 11 of Figure 3, according to the present invention. The higher Q c and Q t values 
of the loudspeaker 20 make the relationship between the reduction in the displacement 

25 response and the increase in co, less straightforward, particularly near the resonance 
frequency co c . To solve this problem, the value of Q c may be "artificially" decreased. 
This is done by setting the value of Q c in Equation 1 to the value of Q c « 6.4 (instead 
of 1/V2 ). Figure 5a shows an example among others of response curves of the low- 
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frequency shelving and notch filter 11 (with a notch at co c by setting Q c =6.4) for an 
under-damped loudspeaker 20, according to the present invention. As can be seen 
from Figure 5a, the resulting response has a notch at the resonance frequency co c , 
which comes from setting the numerator Q- factor in Equation 1 to a value higher than 
5 l/V2 . For this reason, it the filter 1 1 is referred to as the low frequency shelving and 
notch (LFSN) filter. 

The effect of the LFSN filter 11 on the displacement response of the under- 
damped loudspeaker 20 is demonstrated in Figure 5b. The broken line shows the 
loudspeaker's displacement response without the LFSN filter. 

10 The transfer function describing the ratio of the vibration displacement to the 

input signal 22 is a product of the LFSN filter 11 response (transfer function) and the 
loudspeaker 20 displacement response. This is an equalized displacement response 
in the s-domain given by 

H DPE (s) = H c (s)X mVc (s) 

</> 0 s 2 +sa> c /Q c +a> 2 c 1 (2), 

m,R eb s 2 +so) t /Q t +g>? s 2 +sg) c /Q c +g> 2 

1 5 which reduces to 

H DP . E (s) = - A 1 r (3), 

m t R eb S Z +S<£> t /Q t + CD, 

wherein <|)o is a loudspeaker's transduction coefficient (B- 1 factor), Reb is a dc- 
resistance of the voice coil of the loudspeaker 20 and m t is a total moving mass. 

The reduction of Equation 2 to Equation 3 is an important result for operating 
20 the displacement predictor block 14a of Figure 2a. The input to the displacement 
predictor block 14a is the input signal 22, not the output signal 24a from the LFSN 
filter 11 (as in the prior art, see Figure la). Thus the displacement predictor block 14a 
must account for the effect of the LFSN filter 11. It would at first seem that the 
displacement predictor would need to account for the second-order system described 
25 by the loudspeaker displacement response X (s) and the second order LFSN filter 

11, resulting in a fourth-order system altogether. However, the reduction of Equation 
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2 to the single second-order transfer function described by Equation 3 shows that the 
displacement predictor block 14a needs only be a second-order system. 

The same reduction can be made for the z-domain transfer function describing 
a digital processing implementation of the equalized displacement response. The 
product between the z-domain transfer functions of the digital processing version of 
the LFSN filter 11 and a digital model of the loudspeaker 20 displacement is given by 



1 + ^z" 1 + fr 2 . c z" 
" c l + a u z~ l +a 2 , t z~ z \ + a Vc z~ v + a 2 . c z~ 



H DPE {z) = o c o x ^ ; _ 2 1 _ (4), 



wherein <y c is a characteristic sensitivity of the LFSN filter, a x . Vc is a characteristic 

sensitivity of the digital displacement predictor block 14a, b]. c and b2. c are 
10 feedforward coefficients defining the target zero locations, ai , t and a 2 . t are feedback 

coefficients defining the target pole locations and ai- c and a 2 c are feedback coefficients 
defining the loudspeaker's pole locations. 

It is noted that the coefficients bi. c and b2. c can have the same values as ai. c 
and a 2 . c > respectively. Therefore Equation 4 reduces to 

15 H DP . E {z) = a e a — (5). 

l + a u z +a 2 . t z 

The Equation 5 can be written with a single characteristic sensitivity by 
defining 

Vdp_m=°c<*xv c (6), 

wherein ad P _ m is the metrically correct characteristic sensitivity, given by 

R eb k t l -Kc+ b 2c 

wherein ag is a gain of the power amplifier 18a and D/A converter (not shown in 
Figure 2a but used in a case of the digital implementation) and k t is a total stiffness 

of the loudspeaker 20 suspension (loudspeaker's suspension stiffness) including 
acoustic loading from any enclosure. 
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The LFSN filter 11 achieves limiting the vibration displacement by increasing 
the frequency cd,. As shown in Figures 3 and 5 a, increasing this frequency co, reduces 
the gain at lower frequencies, and leaves it unchanged at higher frequencies. This 
provides the desired limiting effect, by changing the displacement response as shown 
in Figures 4a and 5b. 

The displacement-limiting algorithm is shown in more detail in Figure 2b. A 
peak detector 16a-l, in response to the displacement prediction signal 26a from the 
displacement predictor block 14a, provides a peak displacement prediction signal 21 
to a shelving frequency calculator 16a-2. The peak detector provides an absolute 
value of the displacement. It also provides a limited release time (decay rate) for the 
displacement estimate. 

As discussed above, at low frequencies, the gain of the filter varies according 
to the square of the shelving frequency. Due to the nature of the displacement 
response of the loudspeaker 20, it is assumed that the signals that are responsible for 
the excess displacement are at the low frequencies. With this assumption, the 
required shelving frequency is calculated from the excess displacement as follows: 



wherein f r is a shelving frequency required to limit the displacement, f t is a target cut- 
off frequency, xi m and x pn [n] is a displacement predicted by the displacement predictor 
block 14a and normalized to a maximum possible displacement Xmp. 

The maximum possible displacement Xmp can be determined from an analysis 
of the displacement predictor block 20. It can be calculated as 



if{x pn [n] >x lm ) 




(8), 



else 



fr=f, 



mp 



k i R eb 



(8a), 
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wherein gRx is a maximum possible voltage that the D/A and power-amplifier (the 
D/A conversion is used for the digital implementation) can create, and F(Q C ) is a 
function of the loudspeaker's £?- factor, given by 



F(Q e ) = 



1 



i 



(8b). 



10 



15 



20 



The peak value is determined according to 
if ( |xfc["J>*j» [«-!]) 

eke 



(8c), 



wherein x/ w [«] is an instantaneous unity-normalized predicted displacement, jc^f/i] is a 
peak-value of the unity-normalized predicted displacement, and t r is a release time 
constant. The release time constant t r is calculated from the specified release rate d in 
dB/s, according to 

10 - W 



(8d), 



wherein F s is a sample rate. 



The required shelving frequency f r is given by the algorithm of Equation 8. If 
the predicted displacement is above the displacement limit (according to a 
predetermined criterion), this required shelving frequency is increased from the target 
shelving frequency f t according to the first expression of Equation 8. Otherwise (if 
the predicted displacement is below said limit), the required shelving frequency 
remains the target shelving frequency (see Equation 8). If the required shelving 
frequency changes, new values for the coefficients a\. t , a 2 . t , and a c need to be 
calculated by a sensitivity and coefficient calculator 16a-3, thus providing said 
parameter signal 28a to the variable LFSN filter 11. In theory, these parameters could 
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be calculated by formulas for digital filter alignments. However, these methods are 
generally unsuitable for a real-time, fixed-point calculation. Methods for calculating 
these coefficients with polynomial approximations suitable for the fixed-point 
calculation are presented below. 

5 An initial simplification can be made for the f r calculation using Equation 8 by 

defining xi mg , the inverse of the scaled displacement limit, as 

X lmg = l / X lm (9). 

This value, xi mg , is the maximum attenuation needed for the displacement 
limiting. Substituting xim g into the first expression of Equation 8 results in the 
10 following expression for calculating f r : 

fr=f t ^4* P nW (10). 

This value off r is used to calculate co r z , a frequency required for the 
displacement limiting, in rad/s, normalized to sampling rate as follows 

°>,*=^fr (11), 

15 wherein F s is a sampling rate. 

Combining Equations 1 1 and 12 results in 

s 

By defining co,. z in terms off, as in Equations 1 1 and 12 reduces it to 

<>V 2 =yl^l 2 x lmg x pn [n] (13). 
20 From this value of (o r . z , new values of a\. r and a 2 r can be calculated as follows 



(14), 

wherein £ r is a damping ratio. 
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The coefficients ai r and a2 r can be calculated directly from x pn [n], defined as a 
displacement normalized to the maximum possible displacement (x mp ) at a time 
sample n, by combining Equations 10 through 14. Furthermore, these coefficients can 
be approximated by these polynomial series in x pn [n]. 

5 3l r (x p nW)= Pa x 0 + Pa^pn M + P a x 2*\n M + P a x ^\n M + P a^pM (15) 

and 

"l-Axpn [»]) = Pa 2 0 + / Vl*„M + Pa 2 2 x2 pnW + P« r 3^ M + ^ 2 -4^ W (16). 

The characteristic sensitivity a c can be calculated from <5 lr and <i 2 . r according to 

o c =b d (\-a Vr +a 2r ) (17), 
10 wherein 

h *= x h \ h (18). 
l-^,. c +b 2 . c 

The variables 61 c and b 2 c are known from the properties of the loudspeaker 20. 

As b\. c and b 2 . c change only with the loudspeaker 20 characteristics, and thus 
change only infrequently, it is more efficient to compute bj, and store this in a 
1 5 memory for calculating a c . Therefore, according to the present invention, the value of 
bd can to be calculated only once (and not continuously in the real-time), 

The complete formulas for ay r and a 2 r are difficult to approximate with short 
polynomial series for the full range of theoretically valid values of co r . z with an 
adequate accuracy. Potentially, the approximation accuracy can be improved by 
20 increasing the order of the polynomial series. This has not been found to be feasible, 
because it not only increases significantly the complexity of the calculation, it also 
leads to coefficients to be poorly scaled, making them unsuitable for the fixed-point 
calculation. 

The solution to this problem is to optimize the accuracy of the polynomial 
25 coefficients which can mean that different polynomial coefficients will have to be 

used for different hardware and sampling rates, as the latter can be known for a given 
product, so such coefficients can be stored in that product's fixed ROM. 
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Using x p „[n] as the input to the polynomial approximation has an additional 
advantage. Since all ofx pn , a v rl2 7 a 2 - r , and a c are limited to the range (0, 1), the values 
of the polynomial coefficients in the polynomial approximation will be better scaled 
than if, e.g., the required cut-off frequency is used as the input to the polynomial 
approximation Using said x pn [n] simplifies implementation of the polynomial 
approximation using a fixed-point digital signal processor. Therefore, the polynomial 
series can be a good approximation for calculating a\. r and a 2 r from x pn \ 



a l r /2 = - e -W^ cos^T^V^) 



(19), 



a 2 . r = e 



x n „ 

pn 



wherein a f is given by 

and wherein the range of possible values of x pn is 

l) (21X 
This corresponds to a possible range of values of co r z of 

<*rz e(<o t . 2 ,<o t . 2y fi^ ) (22). 

1 5 The Equations 7 through 22 illustrate only a few examples among many other 

possible scenarios for calculating a characteristic sensitivity, ay r and a 2 r by the 
parameter calculator 16a. 

Finally, Figure 6 is a flow chart demonstrating a performance of a signal 
processor with a loudspeaker arrangement utilizing a variable low-frequency shelving 
20 and notch filter 1 1 driven by a feedforward control using a displacement predictor 
block 14a for limiting a vibration displacement of an electro-acoustical transducer 
(loudspeaker) 20, according to the present invention. 

The flow chart of Figure 3 only represents one possible scenario among many 
others. In a method according to the present invention, in a first step 30, the input 
25 electro-acoustical signal 22 is received by the signal processor 10a and provided to 

the LFSN filter 11 of said signal processor 10 and to the displacement predictor block 
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14a of said signal processor 10. In a next step 32, the displacement predictor block 
14a generates the displacement prediction signal 26a and provides said signal 26a to 
the peak detector 16a-l of the parameter calculator 16a of said signal processor 10. In 
a next step 34, the peak displacement prediction signal 23 is generated by the peak 
detector 16a-l and provided to the shelving frequency calculator 16a-2 of said 
parameter calculator 16a. In a next step 36, the shelving frequency signal 23 is 
generated by the shelving frequency calculator 16a-2 and provided to the sensitivity 
and coefficient calculator 16a-3 of the parameter calculator 16a. In a next step 38, 
the parameter signal 28a (e.g., which includes the characteristic sensitivity and 
polynomial coefficients) is generated by the sensitivity and coefficient calculator 16a- 
3 and provided it to the LFSN filter 1 1 . In a next step 40, the output signal 24a is 
generated by the LFSN filter 11. Finally, in a last step 42, the output signal 24a is 
provided to the power amplifier 18 and further to the loudspeaker 20. 

As explained above, the invention provides both a method and corresponding 
equipment consisting of various modules providing the functionality for performing 
the steps of the method. The modules may be implemented as hardware, or may be 
implemented as software or firmware for execution by a processor. In particular, in 
the case of firmware or software, the invention can be provided as a computer 
program product including a computer readable storage structure embodying 
computer program code, i.e., the software or firmware thereon for execution by a 
computer processor (e.g., provided with the displacement predictor block 14a or with 
the parameter calculator 16a or with both the displacement predictor block 14a and 
the parameter calculator 16a). 
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